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Abstract 

A recent high-field magnetization experiment found a phase transition of 
unknown character in the layered, frustrated antiferromagnet RbCuCls, in a 
transverse field (in the layers). Motivated by these results, we have examined 
the magnetic structures predicted by a model of RbCuCls, using the classical 
approximation. At small fields, we obtain the structure already known to 
be optimal, an incommensurate (IC) spiral with wave vector q in the lay- 
ers. At higher fields, we find a staircase of long-period commensurate (C) 
phases (separated initially by the low- field IC phase), then two narrow IC 
phases, then a fourth IC phase (also with intermediate C phases), and finally 
the ferromagnetically aligned phase at the saturation field Hs- The three- 
sublattice C states familiar from the theory of the triangular antiferromagnet 
are never optimal. The C phases and the two intermediate IC phases were 
previously unknown in this context. The magnetization is discontinuous at 
a field « QAHs, in qualitative agreement with experiment, though we find 
much fine structure not reported. 
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I. INTRODUCTION 



The ABX3 family of layered compounds has been studied extensively for several decades 
Much of this interest arises because the frustrated antiferromagnetic interactions in 
the a-b layers give rise to unusual magnetic properties. Quantum and thermal fluctuations 
can have exceptionally large effects in these materials; for example, the stacked triangular 
antiferromagnet (TAFM) CsCuCls exhibits novel fluctuation- induced magnetic phase tran- 
sitions in both a longitudinal field (in the c direction, normal to the layers) and a 
transverse field IQ-Q. 

RbCuCla, another frustrated antiferromagnet of the ABX3 family, is magnetically or- 
dered for temperatures T less than ~ 19K 0; the saturation field Hs (above which the 
ferromagnetically aligned phase is stable) is inconveniently large however (?a 66T) P]. Like 
CsCuCla, RbCuCls is ferromagnetically stacked (in the c direction) and it has an incom- 
mensurate (IC) structure in zero field The Cu^''' ions in the a-b planes of RbCuCls do 
not however form a regular triangular lattice. A room-temperature structural phase transi- 
tion [|l3,|ll|] distorts the structure, yielding a spatially anisotropic intraplane exchange (see 
Figure |I|). As a result, at zero field the spins do not adopt the familiar three-sublattice, 120° 
structure of the TAFM; the structure is instead an IC helical structure along the b direction 
with rotation angle 108° between adjacent spins JT^ . 

A neutron-diffraction experiment |p at T = 1.6K in zero field found that the spins lie in 
the a-b planes and are ferromagnetically aligned in the c direction. The magnetic properties 
in large external field were recently investigated by means of magnetization measurements 
and by ESR A discontinuity was found in the magnetization at a transverse field of 
H = 21. 2T (^ 0.32Hs)- The nature of this field-induced magnetic transition and the 
structure of the high-field state could however not be determined in these experiments. 

The following reports theoretical results for spin structures of RbCuCls in a transverse 
magnetic field. We use the simplest model Hamiltonian and investigate it analytically and 
numerically within the classical approximation. In the field range below Hs, solution of the 
Euler-Lagrange equations yields four IC phases, and in addition many long-period commen- 
surate (C) phases in both the low-field and high-field regions; the three-sublattice C states 
of the TAFM are never optimal. Of these C and IC phases, only the low-field IC (helical) 
phase and the high-field IC (fan) phase were known previously |]13| ; the two intermediate IC 
phases occupy however only a very small field range. The magnetization is discontinuous, 
in qualitative agreement with experiment; the discontinuity occurs however at H = OAlHs, 
rather than at 0.32i^5 as observed, and the magnetization shows additional fine structure. 



II. MODEL HAMILTONIAN 

The model Hamiltonian describing RbCuCls in a transverse magnetic field (which we 
take to lie in the x direction, to be explicit) is ^ 

= — X! ■ Si„+i — ^ 2AJoSfnS-^_^_i — ^ 2JiSin ■ Sj„ — gfisH ^ S'f„. (1) 

in in {ij),n in 

The X and z axes are taken along the a and c directions respectively; Si„ is the spin operator 
(5* = |) of the ith Cu^+ site in the nth a-b plane, and the (ij) sum runs over nearest-neighbor 
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sites in the a-b plane. The first term in Eq. (|l]) is the nearest-neighbor ferromagnetic 
exchange interaction (Jq > 0) in the c direction; the second, with AJq < 0, describes 
an anisotropy of easy-plane type. The third term, with Jj"* < 0, is the antiferromagnetic 
exchange interaction in the a-b planes. As shown in Figure |l|, J^"* is anisotropic; Jl'' = —J[ 
along the b direction while Ji' = — Ji in the other two directions. The fourth and last term is 
the Zeeman energy. The parameters in Eq. (P are Jq = 25. 7K, AJq = — 0.45K, Ji = 10. 6K, 
J[ = 17.4K, and ^ = 2.14 (for H ± c). 

We investigate this Hamiltonian in the classical approximation. The spin operators Sj„ 
then become classical vectors of length i, and the Hamiltonian becomes an energy function; 
the ground-state spin structures are those that minimize the energy. The classical ground 
state of Eq. (|l]) in zero field is the helical structure [ |l^ : 

Sin = 5'X COs(q ■ Yin) + Sj Sm{q ■ Tin) ; (2) 

the easy-plane anisotropy favors spins in the a-b planes. The optimal wave vector is q = 
(0,go)0), with cos go = — Jif^Ji, the above parameter values give go/2vr = 0.2993. 

In the special case J[ = Ji, Eq. (0) reduces to the Hamiltonian of the stacked TAFM. 
Classically, the TAFM ground states are three-sublattice structures; they evolve continuously 
with field H, from the zero- field 120° structure (with go/27r = 1/3 in the above picture) up 
to the aligned state at Hs- These states are continuously degenerate at all H < Hs; energy 
minimization gives only two relations for the three angles (pi between the sublattice spins 
and the field, namely 

cos 01 + cos 02 + cos 03 = qhbH/ {QJiS) (3) 

and I]^=i sin 0j = 0. This degeneracy, which is nontrivial, is broken identically by thermal 
and quantum fluctuations omitted in the classical approximation [0,^. 

We determine the spin structures of RbCuCls in transverse fields with the help of two 
reasonable assumptions: the spins remain in the a-b planes and they are aligned ferromag- 
netically along the c and a directions. We apply periodic boundary conditions in the b 
direction, to reduce finite-size effects. Then, from Figure 1, we can use a chain description 
in which the classical energy per spin is 

1 ^ 

^{{<t^i}) = 7 E f4Ji5' cos(0i - 0,+i) + 2J[S'' cos(0, - (t>i+2) - gfiBHS cos <f)i] ; (4) 
^1=1 

the chain-site index / runs in the b direction, 0; is the angle between S; and the field, and the 
chain length L is a multiple of 3. Equation (^) is just the classical energy per spin for the 
axial next-nearest-neighbour Heisenberg chain, except for the restriction on L. The optimal 
solutions are never TAFM states, but they have related features, and so we use also the 
representation 0j_;/ = 03F+i where i = 1, 2, 3 is the sublattice index. 
The Euler-Lagrange equations are obtained from d£/d(f)i = 0: 

4715"^ [sin(0; - 0i+i) - Sin(0;_i - 0;)] + 2J[S'^ [sin(0; - 0;+2) - sin(0i_2 - (pi)] 

= gfisHS sin 0; . (5) 

These are L coupled, nonlinear difference equations with boundary conditions 0;+^ = (pi 
(mod 2tt). It is not in general legitimate to approximate them by differential equations; 
doing so would miss the sequences of C states that we find at both low and high fields. 
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Equations (j^) have several analytical solutions, namely the zero-field helical state (pi = ql, 
the aligned state 0; = 0, and the TAFM states described by Eq. (||) with Ji replaced by 
J = (2Ji + J'i)/S. Analytical progress is possible for small fields and for H just below 
Hs ||13|. In general Eq. (^) must be solved numerically; this led us to previously unknown 
solutions. The equations have in fact a multitude of solutions. This multiplicity results 
in part from a degeneracy not present in systems (like CsCuCls) where the IC phase is 
driven by a Dzyalshinskii-Moriya interaction. In RbCuCls, both right- and left-handed 
helices are degenerate in energy, at nonzero field as well as at if = (where the relation 
cosgo = —Ji/'^Ji determines only the magnitude of go)- The multiplicity results also from 
the many possible periods for fixed winding number (defined below). 

We solved the equations by Newton's method, starting from approximations obtained 
by various means (for example solutions at a nearby field), using periods L up to ~ 3 x 10^; 
we linearized the equations about the starting values and solved the linear equations to 
obtain the corrections, repeating the procedure until the root-mean-square residual was 
less than 10~^^ or so. The field range from to Hg was covered in increments of AiJ ^ 
2 X l{)~'^Hs- Properties of the IC states (such as the optimal chain length, the energy, and 
the magnetization) are easily obtained by interpolation; for this purpose it is useful to work 
with lengths corresponding to many periods. Periodic boundary conditions necessarily give 
structures with rational periods, but inspection usually suffices to distinguish commensurate 
from incommensurate structures. 



III. MAGNETIC STRUCTURES 

The effect of a transverse magnetic field on helical spin structures arising from competing 
exchange interactions, such as described by the classical Hamiltonian in Eq. (^), was studied 
analytically in detail by Nagamiya et al. They found a distorted helical phase at low 
fields and an IC fan phase at high fields. In the latter, the spins oscillate about the field 
direction; the oscillations decrease with increasing field, vanishing continuously at Hs- If the 
transition from the low-field helical phase occurs directly to the high-field fan phase, then it is 



first-order and occurs at a reduced field estimated as [|T^ Hd Hs = y (1 + P){2 + P) — {1 + (3) 
with P = {1 — Ji/ J[y, or He = 0.4:23Hs for RbCuCls. The structures in intermediate fields 
depend however on the ratio Ji/ J[. For Ji/ J( ^ 1, as in RbCuCla, a TAFM state has lower 
energy than either IC phase for fields near ^Hs- 

Our numerical solution of the Euler-Lagrange equations (^ gives the helical IC and fan 
IC phases at small and large fields respectively, as in Ref. [|T^ . Our results differ however at 
intermediate fields, from H = 0.29Hs to OMHs- For 0.29 < H/Hs < 0.41, we find many 
helical commensurate (C) phases not previously known in this context. These are not TAFM 
states (with period 3 in the chain description) but rather lock-in phases with larger periods; 
they alternate with the helical IC phase, and within our field resolution {AH ^ 2 x lO'^ifg) 
eventually supplant it completely. The behavior is similar in some respects to that found 
in the axial, next-nearest-neighbour Ising model |T^. At H = OAlHs, a strong first-order 



transition occurs to a previously unknown, second IC phase (we call it the IC2 phase); 
this appears to be the transition observed in Ref. [Q. Two weaker first-order transitions 
follow swiftly, to a third IC phase (the IC3 phase, also previously unknown) and then to 
fan phases (both C and IC, intermingled). Finally, there is a second-order transition to the 
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aligned phase; we derive an analytical expression for the saturation field, in terms of the 
parameters in the Hamiltonian. 

Figure our main result, shows the field dependence of the reduced magnetization 

1 ^ 

m = — ^cos0i. (6) 
^ 1=1 

A discontinuity (Am ^ 9 x 10^^) occurs at if = 0A05Hs', this field is however somewhat 
larger than the H = 0.32Hs where a magnetization jump is found experimentally [^, and 
there is a lot of fine structure not reported in Ref. [§. 
We now describe the phases in more detail. 



A. Low-field region: <H < 0.405i?s 

The structures are helical, with three equivalent sublattices; each sublattice phase winds 
through 27r a total of N/3 times, where the total winding number iV is a multiple of 3, like 
the chain length L. 



The helical phase can be treated analytically at small H, as in Ref. 13 . The magnetic 



field distorts the zero-field spin structure 0/ = qqI, modulating the phases. In weak fields, 

ql + a sin ql, (7) 



where the optimal amplitude a and wave number q are field-dependent; they are determined 
so as to minimize the energy S. As the magnetic field is increased, the spin structure 
distorts further as higher harmonics are generated. Figure ^ shows the order parameters 
(the angles (j)i^i> in the three-sublattice representation) found from numerical solution of 
Eq. (§) at ii = 0.33Hs. 

The wave number q of the helical phase is 

' (8) 
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in terms of the total winding number N and the chain length L; one must distinguish between 
q and the wave number 27riV/ L of the chain. /^From Figure |^, at low fields q increases with 
the field, initially quadratically, as the higher harmonics grow. At H = 0.293Hs, however, 
a lock-in transition takes place to a helical C phase with L/N = 13; this phase extends to 
H = 0.296Hs where the helical IC phase becomes optimal again. Other lock-ins occupying 
comparable field intervals occur at L/N = 16, 19, etc in steps of 3, up to L/N = 52 at the 
transition to the IC2 phase. Minor lock-ins (occupying much smaller field ranges) occur at 
L/N = 29/2, 35/2, 18, 20, 41/2, 21, 85/4, 91/4, 23, 47/2, 24, 53/2, 59/2 and 65/2; others 
likely exist but require a finer field step to be detected. Windows of the hehcal IC phase 
exist at higher fields, but they eventually become too narrow to detect (at our field step) 
and we see a staircase of C states which terminates at if = 0.405ii5 in a strong first-order 
transition (Am ^ 9 x 10^"^) to the IC2 phase. 
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B. IC2 phase 



This phase is optimal in the very narrow range 0.405 < H/Hs < 0.419. Figure ^ shows 
the order parameters at if = 0.410ifs'. The structure resembles that of the helical phase 
(as shown in Figure |^); the three sublattices are however no longer equivalent, for the spins 
on only one of them wind through 27r (rather than all three as in the helical phase). The 
optimal chain period increases monotonically with the field, from ^ 180 to ~ 1800 at the 
transition to the ICS phase; the latter transition is weakly first-order, with Am ^ 8 x 10"'^. 



C. ICS phase 

This phase is also optimal over only a very small field range, 0.419 < H/Hs < 0.428. 
The three sublattices are inequivalent, as in the IC2 phase, but here none of the three winds 
through 27r. Figure ^ plots the order parameters at if = 0.424ii5'; one sees that the IC3 
phase is a rippled commensurate phase, in which the order parameter oscillates about a 
particular TAFM state. The ICS state can be expressed in terms of the amplitude and 
the phase 5i of the oscillation on the ith. sublattice as 



+ ttj sin 



2txI' 

+ 5* 



(L/S) 



(9) 



where 0i = 03; the amplitudes are zero for the TAFM state in question. As the field 
increases, the oscillations become smaller and the period L increases (from 80 to ~ 260); 
correspondingly, the TAFM state becomes more favourable, but before it becomes optimal 
there occurs a very weak first-order transition (Am f» 2 x 10~^) to the fan phase. The 
ICS phase (like the IC2 phase) requires a full numerical solution of Eq. (||); it should be 
possible however to demonstrate analytically that the TAFM state becomes unstable to 
these perturbations at some field > 0.428ii5 (the upper limit of the ICS phase). 

We point out a remarkable connection with the theory of fluctuations in the TAFM 



T^,[T5|. Both thermal and quantum fluctuations break the classical degeneracy in the same 
way; for H > ^Hs, two angles are identical in the optimal state, say 0i = 03 0,115 



Our point is that this same state serves as the basis for the ICS state in RbCuCl3; the 
lattice distortion induces the spins to fluctuate spatially, and the energy is minimized when 
the fluctuations occur about the state with two angles identical, just as for thermal and 
quantum fluctuations. 



D. High-field region: 0.428ii5 < H < Hs 

As shown in Figure |^ the spins oscillate about the field direction, without winding 
through 2tt; the three sublattices are again equivalent. The amplitude of the oscillation 
decreases as H increases, going to zero at the transition, at Hs, to the aligned phase. 

The wave number q in this region is obtained from the chain period L as 

g = (27r/S)(l-S/L). (10) 
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As the field increases, the optimal chain length decreases from an estimated value > 10^ at 
the IC3-fan transition, to the interpolated value ~ 29.4 at Hs- Correspondingly q decreases 
monotonically from ^ 27r/3 to ~ 1.88(~ 27r x 0.299); Figure ^ gives numerical values. The 
structure is incommensurate in most of the fan region, but we find also narrow lock-in phases; 
these are less prominent than in the helical region at low fields. The widest C phase, with 
L/3 = 16 and width O.OOGif^, occurs near 0.64if5. Seven others occur at smaller iJ, at 
L/3 = 22, 28 etc; smaller field steps would likely find many more. 

The transition at Hs is second-order and so the fan phase can be treated analytically in 
this region, as in Ref. |]T^. Expanding the energy function 8 m (pi, one finds 

£ = S^AJi + 2J[ - gi2BH/S) 

+^ E {{9f^BH/S - (8Ji + U[)] <P] + + ^J[Ui^2) + 0(0^). (11) 

Introducing the Fourier transform 0i = Z^g^q^*''', one can write (|TT|) as 

E = S'' (4Ji + 2J( - g^^ElS) + {gf^BH/S - 2[J(0) - J(g)]} + 0(0^) (12) 

where J(g) = 4Ji cosg + 2J[ cos2g. Obviously the optimal state of (|12|) is given by 0g = 
for all q if H > Hs = 2[J(0) — J{qo)]S/ g^B, where go gives the minimum value of J{q), 
i.e., cos go = — -^i/2J(. Surprisingly, the optimal IC wave number in the limit H —>■ Hs 
is identical to the value at if = 0, even though the structures are quite different. Our 
analytical result 

Hs = (8Ji + 2 J2/J( + 8J[)S/gfiB (13) 

evaluates to Hs = 66T on using the parameter values given below Eq. (|l|). For H ^ Hs, 
the optimal solution is given by 

(pi ~ acos(go/ + S), (14) 

where 6 is an arbitrary phase; the optimal amplitude a is determined by taking into account 
the fourth-order terms in the energy function S. 

The difference between our result (|13|) and the field (12 Ji + QJ[)S/gfiB (below which the 
aligned state is unstable with respect to TAFM states) is 2 (Ji — J[)'^ S/gfiBJ[- It is then 
obvious that the transition at Hs with decreasing field is to a modulated state, except in 
the spatially isotropic case J[ = Ji. The TAFM states are then rather fragile. They are 
unstable also at if = when J[ ^ Ji, and we show above that they are never optimal for 
the anisotropy appropriate to RbCuCls; it remains to be determined whether they survive 
at intermediate fields for general anisotropy. 

By using the classical approximation, we have omitted the effects of quantum fluctua- 
tions. One can argue that quantum fluctuations are not important for RbCuCls, for their 
effect in CsCuCls in transverse field is to give rise to a plateau in the IC wave number 



T7| , rather than a new phase. On the other hand, quantum and thermal fluctuations in 



CsCuCls stabilize states that are otherwise not optimal, in longitudinal and transverse 
flelds respectively. We made a preliminary effort to understand the effects of quantum 
fluctuations, using a biquadratic term in the energy |T^. Not unexpectedly, we flnd that 



the magnetization develops plateaus, in conflict with experiment, but we did not search 
extensively for other states. 
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IV. SUMMARY 

The classical theory of spin structures of RbCuCla in a transverse magnetic field is 
surprisingly rich. We found four IC phases and as well two series of C phases. We showed 
that the wave number of the helical IC phase at if = is identical to that of the fan IC 
phase as H ^ Hs, and we pointed out a remarkable connection between one of the IC 
phases and the effects of thermal and quantum ffuctuations in the TAFM. We showed that 
the TAFM states of the isotropic case are fragile; infinitesimal anisotropy makes them non- 
optimal at both small and large fields, and they are non-optimal at all fields for RbCuCla. 
We found that the magnetization is discontinuous at ~ OAlHs, in qualitative agreement 
with experiment. Direct observation of the two intermediate IC phases and the C phases 
will be difficult; they occupy only small field intervals, and the large value of Hs requires 
demanding pulsed experiments. Our predictions for the structure of the high-field fan phase 
are likely more easily tested. 
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FIGURES 
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FIG. 3. The order parameter (the angles in the three-sublattice representation, as defined 
below Eq. (4)) in the helical phase at H = 0.33Hs- 
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FIG. 4. The field dependence of the wave number q in the low-field region < H < 0.405Hs. 
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FIG. 5. The angles r in the IC2 phase sX H = QAlHs- 
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FIG. 6. The angles ^j,;/ in in the ICS phase aX H = 0.424ifs. 
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FIG. 7. The angles in the fan phase dX H = 0.66Hs- 
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FIG. 8. The field dependence of the wave number q in the high-field region 0A28Hs < H < Hs- 



13 



